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ABSTRACT 

We present the results of an optical spectroscopic study of 12 GAL£X-discovered star-forming Lya emitting 
galaxies (LAEs) at z 0.3. We measure the emission line fluxes from these galaxies by fitting their observed 
spectra to stellar population models in order to correct for underlying stellar absorption. We revisit earlier 
stellar population model fitting results, finding that excluding now-known AGNs lowers the typical stellar 
population age and stellar mass of this sample to ^ 300 Myr and ^4 x 10^ A^0, respectively. We calculate 
their dust extinction using the Balmer decrement, and find a typical visual attenuation of Ay ^ 0.3 mag, similar 
to that seen in many high-redshift LAEs. Comparing the ratio of Lya/Ha and the Lya equivalent widths to the 
measured dust extinction, we find that the ISMs in these objects appear to be neither enhancing nor seriously 
attenuating the Lya equivalent widths, as would be the case in a quasi-clumpy ISM. Lastly, we perform a 
detailed analysis of the gas-phase metallicities of these galaxies, and we find that most galaxies in our sample 
have Z < 0.4 Zq. We find that at a fixed stellar mass, these low-redshift LAE analogs are offset by ^ 0.6 dex 
lower in metallicity from the general galaxy population at similar redshifts based on the local mass-metallicity 
relationship. This implies that galaxies with Lya in emission may be systematically more metal poor than 
star-forming galaxies at the same stellar mass and redshift, similar to preliminary results at z ~ 2. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 

Galaxies selected on the basis of a bright Lya emis- 
sion line were originally thought to be indicative of primor- 
dial galaxies undergoing t heir first burst of star-formation 
(iPartridge & Peebieslll967b . though they have recently been 
shown to be a complicated group of objects. Studies utiliz- 
ing the technique of spectral energy distribution (SED) fit- 
ting have shown that typical narrowband-selected Lya emit- 
ting galaxies (LAEs) appear to be predominantly young and 
low-mass, with ages < 100 Myr and masses < a few x 
10'' Mr^ (e.g.. iGawiser et al.llMia iFin kelstein et al.J ^MH 



lU" .MfTi (e.g.. lUawiseretalJKUUM l^m Kelstem et al.J IZUU/I 
iPirzkal et al.l2007l:lLai et al.l2007l:lNilsson et al.l2007l) . How 



ever, in many cases these galaxies also appear to contain 
some dust, th us they are hkely not p rimordial in natur e (e.g., 
Pirzkal et al.l l2007t iLai et al I l200a IFinkelstein et"an 120081 
2009d; Pe ntericci et al.l 12009 ). In addition, a small fraction 
of LAEs appear to be more evolved with ages of ^ 0.5 
Gyr and masses of ^ 10'*' Mq, suggesting that there may 
be multiple populations of LAEs, or perhaps a tail in the 
distri bution of LAE properties toward more evolved objects 
(e.g..lFinkelstein etal . 2009d; Lai et al.ll2008l;lPentericci et all 
|2009|). These more evolved LAEs provide a link to the charac- 
teristically more evolved Lyman break galaxies (LBGs; e.g., 
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iKornei et alllMToh . 

While the stellar masses of galaxies can be reasonably 
constrained from SED fitting, the remaining properties suf- 
fer from degeneracies, and can be poorly constrained in the 
absence of rest-frame optical detections, especially at high- 
redshift. Specifically, the derived ages, dust extinctions and 
metallicities are only weakly constrained, as they all result 
in a reddening of the integrated light from a given galaxy. 
While photometry spanning the 4000 A break can improve 
the fidelity of age estimates, the extinction can typically only 
be roughly constrained, and the metallicity not constrained at 
all. In order to obtain more robust estimates of the physi- 
cal make-up of LAEs, direct measurements of the extinction 
and metallicity are necessary. At low redshift, this is typi- 
cally done using measurements of the flux of nebular emis- 
sion lines in the rest-frame optical, e.g., using Balmer line 
ratios to measure the dust extinction, and ra tios of metal lines 
to measure the gas-phase m etallicities (e.g. lKobulnicky et al.l 
ll999tlPettini & Pagelll2004 ). These analyses are not currently 
possible at z > 3, as these diagnostic lines are shifted into the 
mid-infrared. From 2 < z < 3, these measurements are pos- 
sible using near-infrar ed spectroscopy o n 8-lOm cla ss tele- 
scopes (e.g.. lErb et al.ll200 6a: Finkelstei n et al.ll2010l) . How- 
ever, samples of LAEs at these redshifts are only now being 
compiled, and the required integration times are long. 

Presently, we can use low-redshift analogs to try to under- 
stand LAEs at high-redshift. Locally, the Lya properties of 
star-forming galaxies have been studied in great detail. It 
was found that Lya emission can correlate with bright UV 
(and Ha) emission in some regions of a galaxy, but not in 
others. It was also shown that Lya can escape from galax- 
ies even with a significant dust content, especially if outflows 
are present in the interstellar medium (ISM; e.g. Kunth et aT] 
J998i; iHaves et al.ll2"007t lAteFet al. 2008; Ostlin et al. 20090 . 
These results indicate that the mechanisms which regulate 
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Lya escape are complicated, which is an important detail to 
consider when interpr eting the results o f high- redshift LAEs. 

A little further out, iDeharveng et all ( 120081) published the 
discovery of ^ 100 LAEs at z ~ 0.3 using the space-based 
Galaxy Evolution Explorer (GALEX) telescope. A study 
of their luminosity function found that LAEs at low red- 
sh ift are more rare and l ess luminous in L yg than at z > 
3 (IDeharveng et al."2008': 'Co wie et al.ll2010l) . In a followup 
study,rFinkelstein et al. (2009a) studied the stellar populations 
of a subsample of 30 of these LAEs, finding that they ap- 
pear older and more massive than typical high-redshift LAEs. 
A higher fraction of low-redshift LAEs app ear to host active 
galactic nuclei (AGN), from ~ 15 - 40% (iFinkel stein et al.1 
l200Sg iScarlata et alJl2009l: iCowie et alJl2010h . iCowie et al. 
(12010^ compare the metallicities from the ratio of [N II]/Hq; 
emission from z ^ 0.3 LAEs to those of z 0.3 continuum- 
selected galaxies, finding that while the distribu tions over- 
lap, the LAEs extend to lower metallicties. Lastly. lAtek et alJ 
(|2009J used Ha and H/? observations from a subset of these 
galaxies to study the Lya escape fraction, finding clear ev- 
idence for a decrased escape fraction with increased extinc- 
tion, although some galaxies do exhibit a Lya escape fraction 
greater than the continuum, whi ch could imply a clumpy ISM 
geometi-y (e.g. lNeufe ld 1991; Finkelstein et al.ll200 9d). 

In this Paper we present the spectra of our sample of 12 
GALEX LAEs in the Extended Groth Strip (EGS) which 
are dominated by star-formation activity, first analyzed in 
IFinkelstein et al. (2009a) via SED-fitting, and later sh own to 
be devoid of AGN activity in Finkelstein et al.' ('2009b). In §2 
of this Paper, we present the full spectra of every object in 
our sample, along with their measured line fluxes. In §3, we 
revisit the typical ages and m asses of this sam ple, comparing 
to what we earlier derived in IFinkelstein et a l. (2009a) when 
knowledge of AGN activity was unknown. We discuss the 
dust properties and insights into the ISM geometries of our 
sample in §4, and in §5, we present metallicity measurements 
with three separate metallicity indicators, as well as examine 
the mass-metallicity relation, and the implications it has on 
LAEs near and far In §6 we present our conclusions. Where 
applicable, we assume Ho = 70 km s~' Mpc~', il.„, = 0.3 and 
17a = 0.7. 

2. DATA 

2.1. Observations and Data Reduction 

We obtained spectroscopy of 23 of the 27 previously ana- 
lyzed LAEs from Finkelstein et al. ( 2009a) in the EGS, us- 
ing the Hectospec spectrograph with the 6.5m MMT tele- 
scope (all 27 could not be observed in one configuration, 
thus this choice of 23 allowed the most simultaneous ob- 
servations). Hectospec is a multi-fiber spectrograph, with 
300 1.5" fibers covering a 1° d iameter field-of-view at the 
f/5 focus ( Fabricant et alJl2b05h . We acquired spectroscopy 
for our sample using MMT-nHectospec observations, taken 
in queue mode on 19 March 2009, obtaining 120 minutes of 
low-resolution grating spectroscopy, which provides spectral 
coverage from ~ 3650 - 9200 A, with a spectral resolution of 
^ 5 A. The seeing throughout the night was consistently be- 
tween 0.6" - 0.8", and thus much smaller than the Hectospec 
fibers. 

There are two primary pipelines for the reduction of Hec- 
tospec data, IDL and IRAF' based. We chose to use the 

' IRAF is distributed by the National Optical Astronomy Observatory 



IRAF-based pipeline, as it has recently been ported for ex- 
ternal use as External SPECROAD (or E-SPECROAD)^. 
The reduction consists of three tasks: The first, called as 
specroad, makes the master bias, flat and dark frames. The 
next task, specroadcal, traces the fibers across the chips, 
creates a dispersion function, and creates the fiber throughput 
correction. Specroadcal prompts the user to interactively 
fit a wavelength solution, using lines of helium, neon and ar- 
gon. We obtained rms residuals from 0.1 - 0.15 A using 55 
lines spanning the entire wavelength range. Specroadcal 
then applies a dispersion correction to the flat frame. The last 
task, specroadobj, reduces all of the object data, using the 
calibrations made with the previous two tasks, and subtracting 
the sky which was computed from among the six closest of ^ 
50 sky fibers placed throughout the Hectospec field (where 
the flux variation between the fibers was small). The software 
then splits off the 300 fibers into individual fits files, which we 
used moving forward. Each of these files contains the wave- 
length, flux, and flux error array for the specific fiber, where 
the error spectrum is computed from a variance array that was 
created from the raw frame and processed through the same 
reduction steps. 

In addition to our primary targets, we also placed fibers on 
12 F-type stars from the Sloan Digital Sky Survey, which we 
used to flux calibrate the data. These stars are good cali- 
brators, as they are a compromise between strong hydrogen 
lines in hotter stars and other atomic lines in cooler stars. For 
each star, we scaled up a Kurucz (1993) FOV star model to 
match the fluxes, integrating the model spectrum through the 
Sloan Digital Sky Survey (SDSS) g', r' and i' filter curves 
(these three bands were fully encompassed by our wavelength 
range). We then computed a scale factor for each filter by 
dividing the model flux from the SDSS catalog flux for the 
given star A final scale factor was computed as the mean 
of the scale factor from the three filters. This scale was then 
multiplied into the model spectrum. Finally, a calibration ar- 
ray was computed as the scaled up Kurucz model divided by 
the observed stellar spectrum. This was repeated for each star, 
with a final calibration array coming from the mean of each 
of the 12 observed stars. We estimated the accuracy of our 
flux calibration as the standard deviation of the mean of these 
12 observations, which we found to be relatively small, at < 
8%. The final calibration array was multiplied into our object 
spectra, which in addition to flux calibration also corrects for 
the grating transmission function. The full spectra of all 12 
star-forming LAEs are shown in Figure [T| 

2.2. Spectral Fitting 

We measure the emission-line fluxes o f the objects in 
our sample using the method described by 'Moustaka s et alJ 
(2010). B riefly, w e fit a non-negative linear combination of 
Bruzu al & CharlotI ( 120031) population synthesis models with 
instantaneous burst ages ranging from 5 Myr to 10.5 Gyr 
(the age of the Universe at z ^ 0.25), assuming the IChabrieii 
(i2003i) initial mass function from 0.1 - 100 Mq. Motivated 
by the fact that the gas-phase metallicities of these objects 
are sub-solar (see §5), we assume a fixed stellar metallicity 
of Z = 0.004 (2O%Z0); however, we have verified that using 

(NOAO), which is operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation. 

2 E-SPECROAD is available at: 
http://iparrizar.mnstate.ed u/'^juan/research/ESPECROAD/index.php| 
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TABLE 1 

Object Properties and Measured Line Fluxes 



Object 


RA 


Dec 


Redshift Lyo 


[On] 


H/3 


[Olll] 


[Olll] 


Hcj 


[Nil] 




(J2000) 


(J2000) 


A1216 


A3727 


A4861 


A4959 


A5()07 


A6563 


A6583 



EGS7 


215.17155 


53.11394 


0.1996 


16.1 


9.58 ±0.12 


5.29 ± 0.08 


EGS8 


214.43095 


52.46829 


0.2078 


76.7 


2.25 ± 0.08 


0.91 ± 0.06 


EGSIO 


214.55878 


52.39549 


0.2102 


12.3 


4.80 ±0.12 


1.74 ±0.06 


EOS 11 


214.52132 


52.75198 


0.2440 


27.8 


2.05 ± 0.09 


0.84 ± 0.06 


EGS12 


214.70103 


52.29891 


0.2395 


15.2 


2.59 ± 0.08 


1.01 ±0.05 


EOS 13 


214.22613 


52.41111 


0.2462 


35.9 


10.52 ± 0.13 


4.59 ± 0.08 


EGS19 


215.18036 


52.71884 


0.2466 


28.7 


3.13 ±0.08 


1.46 ±0.06 


EGS20 


214.81132 


52.39066 


0.2524 


14.2 


4.65 ± 0.09 


1.53 ±0.07 


EGS21 


215.18599 


52.83513 


0.2515 


25.3 


1.39 ±0.08 


0.77 ± 0.07 


EGS22 


214.30073 


52.59907 


0.2607 


20.8 


1.47 ±0.08 


0.64 ± 0.08 


EGS23 


215.35227 


52.65558 


0.2578 


18.0 


9.75 ±0.11 


4.22 ± 0.09 


EGS25 


214.73180 


52.82437 


0.2630 


20.1 


13.89 ± 0.14 


8.84 ±0.10 



7.79 ± 0.04 
0.86 ± 0.02 
2.01 ± 0.03 
0.28 ± 0.01 
0.70 ± 0.02 
2.73 ± 0.03 
0.47 ± 0.02 

I. 66 ± 0.03 
0.08 ± 0.02 
0.13 ± 0.02 
4.09 ± 0.04 

II. 81 ±0.04 



24.03 ±0.11 
2.52 ± 0.07 
5.93 ±0.08 
1.06 ±0.07 
2.12 ±0.06 
8.43 ±0.10 
1.40 ±0.06 
4.80 ± 0.08 
0.25 ± 0.06 
0.40 ± 0.06 
12.43 ±0.10 
42.91 ±0.17 



19.77 ± 0.20 
2.78 ±0.10 
6.02 ±0.16 
2.98 ±0.10 
3.45 ±0.12 
17.95 ± 0.21 
5.10 ±0.15 
4.43 ±0.11 
2.62 ±0.16 
2.93 ±0.18 
15.27 ± 0.20 
31.02 ± 0.24 



0.89 ±0.11 

< 0.10 

< 0.12 
0.84 ±0.11 
0.44 ± 0.07 
2.59 ±0.17 
1.51 ±0.18 

< 0.21 
1.23 ±0.11 

< 0.14 
1.66 ±0.14 
1.66 ±0.13 



Note. — All line fluxes are in units of 10 " erg s ' cm -. When [N II] is not detected, we report th e 3 cr upper limit, as measured from placement of 
mock lines in the spectra. The Lyo line fluxes are measured from GALEX spectra in Deharveng et al] 1200 8ft (which have a characteristic uncertainty of 
4 X 10"'^ erg s"' cm"-). The redshift is the weighted mean of all lines detected at > 3 cr significance, with a characteristic eiTor Az ~ 0.00001. The last 
column is the best-fiitting extinction on the stellar continuum from the spectral fitting described in §2.2. 



solar-metallicity templates has a negligible effect on the de- 
rived line-strengths. After subtracting the best-fitting stellar 
continuum from the data, we are left with a pure emission- 
line spectrum self-consistently corrected for underlying stel- 
lar absorption. Finally, we fit all the emission lines simulta- 
neously using Gaussian line-profiles. In Figure |2] we illus- 
trate our fitting method for one object in our sample, EGS19. 
In the top panel we plot the observed spectrum (light grey), 
the best-fitting continuum model (black), and the residual 
emission-line spectrum (dark grey). We show the full best- 
fitting emission-line model spectrum (dark red) in the top 
panel, and in the lower panel we zoom into the emission lines 
of interest. The lower, middle panel in particular emphasizes 
the importance of correcting the H/3 line-strength for under- 
lying stellar absorption using this kind of technique (see, e.g., 
iWalcher etal.ll20 1 0, for a recent review). 

The [N II] emission line was not detected in four of our ob- 
jects. In these cases, we computed upper limits to their line 
fluxes by adding in mock emission lines to the spectrum, and 
repeatedly decreasing their line strength until the signal-to- 
noise dropped below 3 a. We found that on average, the 3 
a line flux limits were « 1.0 x 10"'^ erg s"' cm"^. In the 
analysis below, we use the 3 a upper limits as the line fluxes 
for these lines. Additionally, the [O III] A4959 line was not 
detected in two objects (EGS21 and EGS22). This is under- 
standable as these two objects had the lowest [O III] A5007 
fluxes, though both were detected at > 4ct. We thus set the 
[Olll] A4959 flux to be 0.336 x the [Oiii] A5007 flux (the 
theoretical line flux ratio from Storey & Zeippen (2000)) for 
these two objects. We note that out of the remaining 10 ob- 
jects, all but two have their ratio of [Olll] flux ratios at < 3a 
of the theoretical value. The line fluxes for the measured lines 
which are relevant to our analysis are given in Table 1 . 

3. THE AGES AND MASSES OF STAR-FORMING LOW-REDSHIFT 

LAES 

In iFinkelstein et all (I2009a^ we derived the stellar popula- 
tion properties of the parent sample of 30 LAEs via spectral 
energy distribution (SED) fitting, from which we drew the 
23 objects observed with Hectospec, including the 12 con- 
firmed star-formation-dominated galaxies presented here. As 
this was prior to the spectroscopic observations presented in 
this Paper (and .Finkelstein et al.. (.2009b.) ). no knowledge of 



AGNs or line emission was accounted for in the fitting pro- 
cess. The model templates used in this procedure account 
only for stellar emission, thus both AGN continuum emission 
and nebular line emission can cause incorrect fits. Although 
in this paper we exclude the known AGNs from our analysis, 
it is still prudent to revisit the SED-fitting results as the stellar 
masses are key to our analysis below. 

We have re-perfor med the SED fitt i ng usin g the same 
procedure outlined in IFinkelstein et alj (I2009al) . However, 
here we use the u pdated 2007 version of the models of 
'Bruzu al & Chariot ([2003"). In brief, we generate a suite of 
metals, varying the metallicity (0.02 < Z < 2.5 Zq), star- 
formation history (exponentially decaying with 10^ < tsfh < 
10^^), stellar population age (10 Myr - t„„,V(,„f [z]), dust ex- 
tinction (0 < Ay < 1.6 mag) and the ISM dumpiness param- 
eter q (0 < q < 10); where q = TiyalTcontinuum- We fit this suite 
of models to our observations; GALEX far-ultraviolet and 
near-ultraviolet (FUV and NUV, respectively), CFHT Legacy 
Survey u', g', r', i' and z', and the GALEX Lya line flux. We 
find the best-fit model by comparing flux ratios of all bands to 
the u'-band (i.e. colors), and the calculate the stellar mass of 
the best-fit model by a weighted mean of the observed fluxes 
to the fluxes of the best-fit model. Uncertainties on the best- 
fit parameters (excepting stel lar mass) were comput ed via a 
Bayesian likelihood analysis jKauffmann et alj|2003h . where 
the probability distribution is related to the full array via P 

oc e~^ Z^. The array is five-dimensional - one for each of 
our free parameters (5 metallicities x 5 SFHs x 124 ages x 
21 values of extinction x 10 values of q). Using age as an ex- 
ample, the probability of a given object having a given age is 
simply the (normalized) total of P over all other free parame- 
ters. As the mass is determined after the selection of the best- 
fit model, we determine the uncertainty on the masses via 10^ 
Monte Carlo simulations. The emission-line subtracted object 
magnitudes, best-fit stellar population ages, stellar masses and 
dust extinction, as well as the 68% confidence ranges on these 
parameters are given in Table 2. 

Histograms comparing our best-fit stellar population ages 
and stellar masses are show n in Figure [3l cornpared to the re- 
sults for the full sample from lFinkelstein et al.ld2009ah . While 
our earlier study found a number of objects with ages in ex- 
cess of a Gyr, those primarily appear to be dominated by 
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Fig. 1 . — The full Hectospec spectra of the 12 star-fomiing LAEs in our sample. The vertical axis is in units of 10 erg s ' cm ^ A ' , and the reader should 
note that the vertical scale varies from spectrum to spectrum. The typical observed Une width is ~ 6 A. In some cases the 5577 A atmospheric Une is poorly 
subtracted, and should be ignored. 
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Fig. 2. — An example of the emission-line fitting method on EGS19. The 
top panel shows our observed spectrum in gray (smoothed by a 5-pixel box- 
car), with the best-fitting spectral model shown as the black curve. The un- 
derlying emission line spectrum, after subtracting the stellar continuum, is 
shown in red. The bottom three panels show the regions around some of the 
lines of interest to our study. Of note is the strong stellar Balmer-line ab- 
sorption. Connecting for this underlying absorption is crucial to an accurate 
measurement of the Balmer decrement. 

AGN-hosting LAEs, thus the ages may be unreliable. We find 
that the updated ages for our star-forming LAE sample are 
primarily between 100 - 1000 Myr, with a median of ^ 300 
Myr We find a similar result for the stellar masses, shown in 
the right-hand plot of Figure [3] in that the many LAEs with 
derived stellar masses > lO"' Mq appear to be dominated 
by AGN hosts, and that the star-forming LAEs in our cur- 
rent sample have masses ^ 10'* - 10'° Mq (with a median 
of 4 X 10^ Mq). While these low-redshift LAEs are older 
and more massive than typical narrowband selected LAEs, 
they are similar to those of the more evolved, rarer, LAEs 
seen at hig h redshift (e.g.. [ Finkelstein et al.ll2009dt iLai et al.l 
HOOS; Pent ericci et al.ll2009h . However, we do caution that 
the derived ages are highly uncertain, as although many ob- 
jects have best-fit ages on the order of a few hundred Myr, the 
uncertainties on these derived ages are high - sometimes or- 
ders of magnitude. As shown in Table 2, the uncertainties on 
the stellar masses are much smaller, and we use these best-fit 
masses with their uncertainties in our metallicity analysis in 
§5. 

4. DUST EXTINCTION AND THE ISM 

Using the Balmer decrement we can measure the level of 
dust extinction in our galaxies, providing a much more robust 
estimate than that derived from SED fitting (listed in Table 
2), which is subject to a number of degeneracies. Tradition- 
ally, this is done by measuring the ratio of Ha/H/3, and then 
comp aring to the theo retical ratio, which we take as Ha/H/3 = 
2.86 ( Osterbrockl 1 1 989i) . Any increase in this ratio over the 
theoretical value is likely due to dust, and, assuming an ex- 
tinction law, one can compute the color excess. Table 3 tabu- 
lates the extinction results assuming the Calzetti et al. (2000) 
starburst extinction law. The computed color excesses for our 
objects range from < E(B-V)^ < 0.4, with 10/12 objects 
having E(B-V)j significantly greater than zero (at la signif- 
icance; 9/12 at 2a). We converted the color excess on the 
nebular emission to an extinction on the stellar continuum as- 



suming that E(B-W)steiiar = -44 E(B-V)i,,„, as has been found 
locally dCalzetti et al.ll2000l) . We then converted to the stellar 
attenuation in the V-band via the adopted attenuation curve. 
The derived attenuation on the stellar continuum ranges from 
Ay = - 0.7, with a mean Ay ^ 0.3. Comparing the de- 
rived attenuations from the Balmer decrement to those from 
the photometric SED fitting (listed in Table 2), we find excel- 
lent agreement, as the difference in attenuation for 9/12 galax- 
ies are smaller than their combined la uncertainties. This 
implies that in this sample of galaxies the conversion used 
from gas extinction to stellar extinction is accurate, and thus 
that the nebular emission lines are more extincted than the 
stellar continuum. However, the uncertainties on the Balmer- 
decrement-derived attenuations are much smaller than those 
fr om SED fi tting, highlighting the advantage of spectroscopy. 

lAtek et al . (2009) studied the Lya escape fraction in LAEs 
from the z 0.3 sample of IDeharveng et al.l (1200 8). obtain- 
ing spectroscopic measurements of objects from the Chandra 
Deep Field - South (CDFS) and ELAIS-Sl fields. As part of 
their analysis, they also computed the dust extinction using 
the ratio of Ha/H/3. From their Figure 2, they found values of 
the color excess from -0.3 < E(B-V) < 0.5, covering roughly 
the same range of color excess valu es as we find in our in- 
dependent sample. lAtek et al.l (l2009h fi nd that 12/21 of their 
objects are consistent with E(B-V) > 0. IScarlata etall (l2009h 
also studied z ~ 0.3 LAEs, obtaining optical spectro scopy 
of 31 galaxies drawn from the iDeharveng et aTl ( 120081) EGS, 
NGPDWS and SIRTFFL samples, with their EGS sample in- 
cluding four objects in our sample of 12 star-forming galax- 
ies. They do not compute the dust extinction, but to facilitate 
comparison we have computed the color excess based on their 
tabulated Ha and H/3 fluxes for the 20 objects which have 
3 a detections of both emission lines. Their objects have a 
range of -0.05 < E(B-V) < 0.74, with 11/20 objects consis- 
tent with E( B-V) > at > 1 0-. While both Atek et al. (2009!) 
and Scarla ta et al.l ((2009) find a slightly smaller fraction of 
objects with significant dust attenuation than we do, this dif- 
ference is likely due to the small number of objects analyzed 
in all three studies, as well as the large uncertainties on the 
extinction in general. We note that Atek et al. ( 2009) used 
the extinction curve of Cardell i et al.l(ll98 9). Had we adopted 
this curve, we would have found higher values of the color ex- 
cess, though the derived attenuation would be similar to what 
is reported in Table 2. 

While many high-redshift LAEs typically have derived 
non-zero extinctions, these are inferred from SED fitting. As 
discussed in Finkelstein et al. (2009a) and §3, the galaxies in 
our sample are more evolved than typical narrowband selected 
LAEs at high redshift. Thus, it is interesting to compare the 
mean value o f Ay ^ 0.3 mag mea sured here to high-redshift 
LAEs. While iGawiser et an(l2006h derived no dust extinction 
in LAEs at z ^ 3.1, others have frequently found evidence 
for du st in typical narrowband-selected LAEs. Pirzkal et aD 
(120071) found Ay = - 0.6 mag in LAEs at z 5, Lai et all 
(20071 found Ay = 0.4 - 1.7 mag in LAEs at z ~ 5.7, and 
fFinkelstein et al.1 (l2009dh found Ay = 0.1 - 1.1 mag in LAEs 
at z ~ 4.5. Thus, though it appears that although high-redshift 
LAEs are on average younger and lower mass, many of them 
appear to contain a similar range of dust attenuation as these 
low-redshift LAEs. 

4.1. ISM Geometry 

The geometry of the ISM can play a strong role in the es- 
cape of Lya photons. Lya photons are resonantly scattered 
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Fig. 3. — Left: Distribution of best-fit stellar populations ages for the 27 LAEs from 'Finkelstein et al.' 12009a) are shown as the gray histogram. The best-fit 
ages we derive here for the subsample of 12 star-forming LAEs are shown as the blue histogram. While the original sample of 27 LAEs had many objects with 
best-fit ages > 1 Gyr, when we exclude now-known AGNs from Finkelstein et al. ( 2009b), we find that all but one galaxy appears younger than 1 Gyr The dotted 
blue histogram denotes the ages for the same 12 LAEs from Finkelstein et al. (2009a), which broadly agree with our updated results. Right: Same as the left, for 
the derived stellar masses. Similar to the ages, a smaller characteristic mass appears likely for these low-redshift LAEs once known AGNs are excluded. 



TABLE 2 
SED Fitting Results 



Object vafuv 


mNuv 


mi, mg m,- m, mj f'Lya/lO '"^ 


Age 


Age 


Mass 


Mass 


Av 


Av 


(mag) 


(mag) 


(mag) (mag) (mag) (mag) (mag) (erg s"' cm^^-. 


Best-Fit 


68% Range 


Best-Fit 


68% Range 


Best-Fit 


68% Range 








(Myr) 


(Myr) 


(10' Al0) 


(IC Mq) 


(mag) 


(mag) 



EGS7 


21.63 


21.29 


21.31 


20.88 


20.57 


20.45 


20.48 


2.8 


450 


260- 


800 


1.3 


1.2- 


1.5 


0.32 


0.08 - 


0.57 


EGS8 


22.03 


21.66 


21.88 


21.18 


20.76 


20.64 


20.63 


1.5 


320 


140- 


400 


1.6 


1.6- 


2.1 


0.00 


0.00- 


0.40 


EOS 10 


21.59 


21.29 


21.45 


21.02 


20.75 


20.64 


20.60 


2.9 


570 


320- 


900 


1.6 


1.3- 


1.7 


0.32 


0.00- 


0.49 


EOS 11 


22.16 


21.77 


21.34 


20.58 


20.07 


19.85 


19.73 


1.4 


290 


80- 


360 


3.8 


3.7- 


4.9 


0.49 


0.32- 


1.05 


EOS 12 


21.64 


21.35 


21.14 


20.42 


19.93 


19.76 


19.64 


2.5 


320 


90- 


2600 


4.9 


4.8- 


5.8 


0.49 


0.16- 


0.89 


EOS 13 


21.76 


21.23 


21.24 


20.61 


20.13 


19.93 


19.80 


2.1 


810 


400- 


2600 


4.2 


3.8- 


7.3 


0.57 


0.16- 


0.81 


EOS 19 


21.86 


21.29 


21.28 


20.53 


20.03 


19.82 


19.67 


2.2 


320 


180- 


4200 


4.9 


4.8- 


10.8 


0.57 


0.16- 


0.89 


EGS20 


22.02 


21.64 


21.71 


21.21 


20.81 


20.74 


20.66 


1.8 


720 


230- 


1600 


2.0 


1.7- 


2.0 


0.32 


0.08- 


0.49 


EGS21 


21.75 


21.28 


21.11 


20.30 


19.58 


19.29 


19.11 


1.9 


9000 


2600- 


8800 


34.4 


10.5- 


33.7 


0.00 


0.00- 


0.65 


EGS22 


21.99 


21.49 


21.32 


20.62 


20.02 


19.82 


19.67 


1.7 


260 


160- 


3500 


4.8 


4.7- 


4.9 


0.65 


0.24- 


0.89 


EGS23 


21.88 


21.37 


21.50 


21.15 


20.82 


20.63 


20.55 


2.9 


80 


27- 


570 


1.1 


1.1 - 


2.0 


0.97 


0.32- 


1.38 


EGS25 


22.18 


21.10 


21.44 


21.17 


20.95 


20.87 


20.65 


2.6 


60 


18- 


71 


1.2 


1.2- 


2.2 


0.97 


0.89- 


1.46 



Note. — Th ese magnitudes have had all emission lines detected at > 3a- subtracted. The magnitude errors are 0.2 and 0.02 mag in the GALEX and optical bands, 
respectively (see lFinkelstein et"al]42009"3) ). 



by H I, thus in a homogeneous geometry (i.e. a dust screen in 
front of the sources), where H I gas and dust are evenly mixed, 
Lya has a difficult time escaping, and thus has a much greater 
chance of encountering a dust grain than the continuum pho- 
tons. On the other hand, if the ISM is inhomogeneous, with 
the dust locked together with neutral gas in many clumps re- 
siding in a tenuous, ionized medium, the Lya photons will 
resonantly scatter at the edge of these clumps and avoid 
the dust, while the continuum light will pass throu gh and 
be subject to dust attenuation (Neufeld 1991; Hansen & OhI 
[2006). This has the net effect of increasing the Lya escape 
fraction over that of the continuum, and thus increasing the 
Lya equivalent width (EW). This has been used to explain 
the large number of high Lya EW LAEs which have been 
fo und (e.g., Kudritzki et al. 2000; Malhotra & Rhoads 2002). 
In lFinkelstein et all (i2008i) and iFinkelstein et al.i (2009d), we 
included this scenario in our stellar population models. We 
parameterized the ISM geometry with a free parameter "q", 
which is defined by 

flya = flya X exp(-^rc) ( 1 ) 



where is the optical depth of dust on the continuum light, 
thus q = Tiya/Tc- We make the assumption that the cov- 
ering fraction of clumps is large, meaning that a Lya pho- 
ton will scatter many times prior to escape. In this sce- 
nario^ the extinction curve (which we assume to be that of 
ICalzetti et al. (2000)) is the same as that for a uniform dust 
screen (Scarlata et al.ll200"9h . excepting that Lya photons can 
preferentially avoid encountering the dust. When ^ = 10, this 
has the effect of a homogeneous ISM, where Lya is attenu- 
ated far more than the continuum as resonant scattering in- 
creases the chance of absorption by dust. When q = Q, this 
simulates the Neufeld (1991) scenario, where all of the dust 
is locked away inside neutral gas clumps. In between values 
of q reflect intermediate geometries, with any value of <7 < 1 
enhancing the Lya EW. Using this, we found that dust en- 
hancement can explain the SEDs of many z ~ 4.5 LAEs, as 
10/14 objects had q <\. We also used the ISM du mpiness 
parameter q when we fit models in IFinkelstein et al.l (t2009ah 
and §3 to the present sample. Our 12 low-redshift LAEs have 
a best-fit range from < ^ < 1, with 5/12 objects having q 
from 0.75 -1.0 (using the updated model fitting results dis- 
cussed in §3), reflecting an ISM which does not much enhance 
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or attenuate the Lya EWs. 

With these spectroscopic data, we can now perform a more 
in depth analysis into the ISM geometry by comparing the 
resonantly scattered Lya line to a non-resonant line, such as 
Ha, as a proxy to trace the continuum photons. Under Case 
B recombination, the intrinsic ratio of Lya/Ha is ^ 8.7. Any 
deviation from this value is presumably due to dust (although 
preferential scattering of Lya photons could play a role). In 
a typical homogeneous ISM, this value will drop quickly if 
dust is present. If, however, the ISM is in a configuration 
that results in the enhancement of the Lya EW, this value will 
increase. 

Our Ha measurements come from a 1.5" diameter fiber, 
while the Lya measurements are from the far-ultraviolet 
(FUV) GALEX prism, which h as a beam full-w idth at half- 
maximum (FWHM) of 4.2" (iMorrissey et al.l ^007). As 
many of the se galaxies are larger than 1.5" (see Figure 1 of 
iFinkelstein et al. (2009a)), we need to account for aperture 
differences between the Lya and Ha measurements. We did 
this by first measuring photometry of our sample of galax- 
ies in the Canada-France-Hawaii-Telescope Legacy Survey 
(CFHTLS) i' image (which has a seeing of ^ 0.7", compara- 
ble to our Hectospec observations), using an aperture diame- 
ter of the same size as the Hectospec fibers (L5"). We then 
smoothed this image such that the FWHM of stars in the im- 
age matched the FWHM of the GALEX beam, and performed 
photometry on our sample again, this time using 4.2" diam- 
eter apertures. In both cases, the photometry was done us- 
ing the Source Extractor software package (Bertin & ArnoutsI 
[1996) . The aperture correction which we applied to the Ha 
flux measurement is thus the ratio of the flux in the 4.2" aper- 
ture to that in the L5" aperture. This aperture correction fac- 
tor ranged from 1-4, with a median of 1 .3, and 1 1/12 galaxies 
having an aperture correction factor of < 2.0 (one object had 
a correction of < 1, which we rounded up to 1). 

The top panel in Figure|4]plots the Lya/Ha ratios, using the 
aperture correction, for our sample versus their derived dust 
attenuation from the previous section. Examining first the 
Lya/Ha ratios, we find they span 1 < Lya/Ha < 7, with all 
12 objects having line ratios significantly less than 8.7, imply- 
ing dust is likely present. However, the exact amount of dust 
cannot be derived from this ratio alone as the geo metry plays 
a role. U sing stellar population models from Fink elstein et al.l 
(l2009al) . we also plot tracks of constant ISM dumpiness, from 
q = - 10. Thus, by examining where our objects fall in this 
plane, we can gain insight into their ISM geometry. 

Investigating this plane, we find that the majority of our ob- 
jects are consistent with values of q ~ 1-3. This implies that 
the dust that is present lies in a geometry somewhere between 
that of a uniform screen and that of purely clumps. In the 
bottom panel of Figure |4] we plot the GALEX derived Lya 
EW versus the attenuation we derive here. We see that the 
Lya is between 20 - 40 A for most objects regardless of the 
level of extinction. This again implies that the ISM is neither 
overtly enhancing nor depleting the Lya EW (similar also to 
the resul ts from SEP fittin g). Similar results are seen at z '--^ 
2, where iBlanc et al.l (1201 0) recently analyzed a sample of 
100 LAEs and found that Lya and the rest-frame UV contin- 
uum experienced similar levels of extinction (i.e. q = 1). As 
mentioned above, Finkelstein et al. (I2009dl) found than many 
z ^ 4.5 LAEs were consistent with q < 1, thus at much higher 
redshift the ISMs may be more clumpy, though this was indi- 
rectly inferred via SED fitting. 

,Atek et al.. (.2009.) studied the escape fraction of Lya pho- 



tons in their sample of z ^ 0.3 LAEs, finding that some objects 
indicate Lya-assisted escape due to a clumpy ISM (q < 0), 
some objects appear to have similar Lya and UV-continuum 
escape fractions (q ~ 1), and some objects have Lya es- 
cape fractions < the UV continuum (q > 1). iScarlata et aP 
(2009) recently investi gated the ISMs of LAEs from the 
IDeharveng et al.l (120081) sample, examining a number of ra- 
diative transfer scenarios, including a clumpy ISM. In partic- 
ular, they study how the hydrogen line ratios change with the 
number of clumps, deducing that they have clumpy interstel- 
lar media. However, these authors also conclude that there 
is no evidence for preferential shielding of Lya from dust, as 
they do not have any objects with Lya/Ha > 8.7, similar to 
our results. Scarlata et al. note that the attenuation curve in a 
clumpy ISM can differ from that for the uniform dust screen 
if the number of clumps are small. Thus although we have 
assumed a large number of clumps such that the uniform dust 
screen law applies in our above discussion, this need not be 
the case. 

Examining the results for all studied GALEX LAEs, it ap- 
pears that although the ISMs in low-redshift analog LAEs 
span a wide range of extinction and dumpiness values, some 
level of clumping appears likely (i.e., a uniform dust screen 
does not appear to be supported by the data). We note that the 
necessity of an aperture correction based on photometric data 
adds a level of overall uncertainty, as it assumes the Ha flux 
is spatially correlated with the broadband emission. This can 
be avoided in future studies when both lines are measured in 
similar apertures. 

5. MASS-METALLICITY RELATION 

Using our measured emission lines, we can measure the 
gas-phase metallicities of the galaxies in our sample with a 
number of metallicity indices. We can then examine these 
galaxies on a mass-metallicity relation to see how they com- 
pare to other populations. We use ratios of emission lines 
of oxygen and nitrogen to those of hydrogen to measure the 
gas-phase metallicities of the galaxies in our sample. Specif- 
ically, we use the calibrations of the N2 and 03N2 indices of 
,Pettini & Pagel (2004) , where 



N2 = log 



/[Nil] 
V Ha 



and 



03N2 = logr[251^^^V 
^ V [Nil] /Ha J 



(2) 



(3) 



Additionally, we obtain another measure of the metallicity us- 
ing the well-known R23 index, where 



^23 ■ 



[Oil] A3727 + [QUI] A4959 , 5007 
up ' 



(4) 



using the calibration of iKobulnickv et al] (119991) . This in- 
dex is double valued above Z ~ 0.3 Zq, thus we tabulated 
the metallicities for both the upper and lower branches. Ta- 
ble 4 presents the measured metallicities of our galaxies from 
each of these indices, as well as their errors, including sys- 
tematic uncertainties from the calibrations of 0.05, 0.14 and 
0.18 dex for R23, 03N2 and N2 indices, respectively. We 
note that these indices were calibrated using nearby galax- 
ies, and they may be sensitive to a change in Hll region 
physical conditions such as the ionization parameter or elec- 
tron temperature. In fact, many high-redshift star-forming 
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Fig. 4. — Left: The distribution of Lya/Ho ratios versus the measured dust attenuation for the 12 star-forming LAEs in our sample. The values of Lyo/Ha 
have been con'ected for aperture differences. The colored Hnes show model tracks for constant values of q, the ISM dumpiness factor By studying our objects 
in this plane, we see that few appear to have an ISM geometry that preferentially favors Lya escape over continuum photons. Most of the sample lies at q ~ 1-3, 
implying that Lya encounters a similar level of dust (and possibly a bit more) than than the continuum. Right: The Lya EW versus the derived dust attenuation. 
The EW is roughly invariant to the level of dust, implying that the ISMs of these objects are neither strongly enhancing nor depleting the Lya EWs, and thus that 
neither a uniform dust screen, or a purely clumpy ISM geometry apply to these objects. 



TABLE 3 

Measured Dust Extinction and Gas-Phase Metallicities 



Object E(B-V),,„, 








z' 


7" 


Z Weighted Mean 




(mag) 


(Z©) 


(Z©) 


(Z©) 


(Z©) 


(Z©) 



EGS7 





23 


±0.02 


0.41 ± 0.04 


0.30 ± 


13 


0.27 ± 


09 


0.18 


± 


0.03 


0.68 


±0 


11 





20 


± 


0.03 


EGS8 





06 


± 0.06 


0.11 ±0.11 


0.26t 




0.29t 




0.20 


± 


0.12 


0.86 


±0 


36 





25 


± 


0.07 


EGSIO 





16 


±0.04 


0.29 ± 0.07 


0.18t 




0.23t 




0.28 


± 


0.11 


0.66 


±0 


20 





22 


± 


0.05 


EGSll 





19 


± 0.07 


0.33 ±0.12 


0.85 ± 


38 


0.73 ± 


26 


0.23 


± 


0.17 


1.11 


±0 


46 





44 


± 


0.13 


EGS12 





15 


±0.05 


0.27 ±0.10 


0.54 ± 


24 


0.48 ± 


17 


0.24 


± 


0.13 


0.90 


±0 


33 





36 


± 


0.10 


EGS13 





27 


±0.02 


0.47 ± 0.04 


0.58 ± 


24 


0.53 ± 


17 


0.24 


± 


0.05 


0.96 


±0 


15 





28 


± 


0.05 


EGS19 





17 


± 0.05 


0.30 ± 0.08 


0.87 ± 


38 


0.81 ±0 


28 


0.17 


± 


0.09 


1.31 


±0 


34 





26 


± 


0.08 


EGS20 





01 


±0.04 


0.02 ± 0.08 


0.3lt 




0.3 It 




0.25 


± 


0.11 


0.73 


±0 


24 





28 


± 


0.07 


EGS21 





14 


±0.09 


0.26 ±0.16 


1.13±0 


53 


1.33 ±0 


52 


0.14 


± 


0.14 


1.74 


±0 


66 


1 


36 


± 


0.33 


EGS22 





40 


±0.11 


0.72 ± 0.21 


0.3 It 




0.52t 




0.38 


± 


0.59 


1.07 


±0 


85 





38 


± 


0.13 


EGS23 





20 


±0.02 


0.36 ± 0.04 


0.49 ± 


21 


0.41 ± 


13 


0.22 


± 


0.05 


0.80 


±0 


14 





26 


± 


0.05 


EGS25 





17 


±0.01 


0.31 ±0.03 


0.33 ±0 


14 


0.28 ± 


09 


0.14 


± 


0.02 


0.73 


±0 


10 





15 


± 


0.02 



Note. — t [n hj was undetected in EGS8, EGSIO, EGS20 and EGS22, and the 3 cr upper limit was used on the Une 
flux to compute N2 and 03N2. N2 is thus an upper limit and 03N2 is a lower limit. The weighted mean metallicity was 
computed with just one value for R23 per object, using the value which gave the lowest error on the weighted mean with the 
N2 and 03N2 metallicities (which ended up being the lower branch in all cases). The exception is EGS21, where although 
the weighted mean uncertainty was formally smaller with the lower R23 result, the high metallicities with the 03N2 and 
N2 indices imply that the upper branch is correct. 



galaxies may have an increased ionization parameter as they 
appear enhanced in [Olll]/H/3 from the l ocal star-forming 
galaxies sequences from the SPSS (e.g.. lErb et alj l2006bl : 
iFinkelstein e t al."2009c': 'HainUne et al]|2009ir However, as 
shown in|Finkelstein et al. (2009b) the star-forming members 
of our sample have similar [Olll]/H/3 ratios as the SDSS 
galaxies, thus they likely have similar H II region properties 
as well. 

Prior to this computation, the emission lines were corrected 
for the inferred dust extinction from §4. We also present a 
weighted mean metallicity from all three measures, where 
the R23 value which results in smaller uncertainty on the 
weighted mean was chosen (with the exception of EGS21, 
where it is obvious that the upper branch of R23 is the ap- 
propriate one). In nearly all cases, the three measures of the 
metallicities are consistent within their uncertainties, thus our 
results are not dependent on the precise indices used. Investi- 



gating Table 3, we find that these low-redshift LAEs typically 
have metallicities from 0.2 - 0.4 Xq. 

In Figure |5] we plot the 12 star-forming z ^ 0.3 LAEs 
in a mass-metallicity plane. Star- forming galaxi e s wer e 
extensively studied in this plane by iTremonti et al.l (120041) . 
who used SDSS spectra of ^ 53,000 star-forming galax- 
ies to study the mass-metallicity relation at z ^ 0.1. In 
that work, they computed the metallicities by using a com- 
bination of all strong emission lines. Because different 
metallicity indicator s suffer inherent systematic uncertain- 
ties o f 0.2 - 0.6 dex dKewlev & Dopitall200l ITremonti et alj 
2004; Moustakas et al.l2010l) . we recomputed metallicities for 
the SDSS sample using the same techniques as above. In 
Figure E] we sh ow the best-fit mass-metallicity relation of 
ITremonti et alj ([2004), which is shifted downward by 0. 1 dex 
to match our recomputed SDSS metallicities. 

At any given mass, all but the most massive of our z ^ 
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Fig. 5. — Our sa mple of 12 star-formin g low-redshift LAE analogs plotted in a mass-metallicity plane. For comparison, we show the sample of ~ 53,000 
SDSS galaxies from'Tremonti et al." (2004"), where we recomputed their metallicities using the same methods as our sample. The solid line is the best-fit mass- 
metallicity relation from Tremonti et al. 1 2004) (shifted down 0. 1 dex to match the sequence shown using our metallicity derivation). The dashed line is the SDSS 
relation shifted downward by ~ 0.6 dex, highlighting that these LAE analogs are significantly shifted downward in metallicity for a given mass when compared 
to the SDSS sample. This implies that galaxies exhibiting Lyo in emission may be systematically more metal poor than other galaxies at the same stellar mass 
and redshift. 



0.3 LAEs lie significantly below the z ~ 0.1 mass-metallicity 
relation. In Figure |5] we show the best-fi t curve t o the z 
'--^ 0. 1 mass-metallicity relation from Tremonti et alJ (l2004h 
shifted to match our derived SDSS metallicities, and we find 
that shifting it downward by a further ~ 0.6 dex best matches 
our sample of LAEs. This offset is larger than plausible er- 
rors in the stellar masses. This illustrates that at low-redshift, 
galaxies exhibiting Lya in emission have significantly lower 
metallicities tha n comp a rable m ass galaxies with no such pre- 
selection. Cowie et al.l (120101) found similar results, as they 
noted that GALEX-selected LAEs had lower metallicities than 
star-forming galaxies selected via the near-UV continuum. 
However, their brighter LAEs had similar metallicities to their 
continuum-selected sample, while all but one of our more 
massive LAEs are still more metal-poor than the SDSS sam- 
ple. Th e differences may b e due to the use of solely the N2 
index in lCowie et al.l ( 1201 0). as it typically yields the highest 
metallicities of the three indicators we use in our sample. 

At z 0.3 the mass-metallicity relation of LAEs is shifted 
lower than that of the SDSS galaxies, which represent the gen- 
eral star-forming galaxy population at these redshifts. At high 
redshift, LBGs ap pear to be the dominant galaxy population. 
lErb et alJ (l2006bh has studied the mass-metallicity relation 
of LBGs at z ~ 2.3, finding that their objects fit the shape 
of the local mass-metallicity relation, but shifted downward 
in metallicity by 0.56 dex (over a metallicity range of 0.4 - 
0.9 Zq for galaxies with stellar masses of 10^ - 10" Mq)- 
Current results show that high-redshift LAEs are on average 
less evolved in age, mass and dust content than LBGs (e.g,. 



Gawiser et aLll2006l; IPinkelstein et"aI1l2009dt iPentericci etaTI 



2009), thus it may seem intuitive that they would be lower 
in metallicity as well (i.e., less evolved in all physical char- 
acteristics), possibly shifted even further down for the LBG 
mass-metallicity relation. 

This was recent ly studied at high-redshift by 
'Finke lstein et al.l (120101) . where upper limits on the metal- 
licities of two LAEs at z = 2.3 and 2.5 were measured to 
be Z < 0.17 and 0.28 Zq, respectively, with the N2 index. 
The object with the lower metallicity upper limit appears to 
be sign ificantly below even the z = 2.3 LBGs of lErb et"al] 
(2006b) on the mass-metallicity relation, analogous to our 
results at z ^ 0.3. At z = 2.3, LBGs represent the general 
star-forming galaxy population, and LAEs appear to be more 
metal-poor than this population (though with a sample of 
two objects we cannot draw robust conclusions). However, 
we see something very similar at z 0.3, where the SDSS 
sample represents the general star-forming galaxy population. 
Thus, it may be that galaxies exhibiting Lya in emission have 
preferentially lower metallicities. However, the evidence for 
this at high redshift is not strong, and a confident resolution 
will not be possible until we can measure the metallicities 
of large samples of high-redshift LAEs, which will be 
done in the near future with the upcoming multi-object NIR 
spectrographs. 

6. CONCLUSIONS 

We present results from a spectroscopic survey of GALEX- 
discovered LAEs at z 0.3. These objects are useful as 
nearby analogs to high-redshift LAEs, as we are able to di- 
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rectly measure physical quantities which are only inferred at 
high-redshift. In this Paper, we present the full optical spectra 
of 12 star-forming LAEs, including hydrogen, oxygen and ni- 
trogen line fluxes. We first re-examine their stellar population 
properties, and find that when excluding the recently identi- 
fied AGNs in the parent sample, the average ages and masses 
of the remaining star-forming LAEs are lower, at ^ 300 Myr 
and 4 x 10"' Mq respectively. These are older and more 
massive than a typical high-redshift LAE, but similar to a rare 
population of evolved high-redshift LAEs. 

Correcting our measured emission line fluxes for stellar ab- 
sorption, we use the Balmer decrement to investigate the level 
of dust attenuation in these galaxies. We find that these galax- 
ies all contain some level of dust, with an average extinction 
of Ay ~ 0.3 mag. This is similar to what is typicaUy found 
at high redshift, and is also consistent with the overall trend 
of the existence of dust in galaxies which are still able to emit 
in Lya. Using the ratio of Lya/Ha and the measured extinc- 
tion as a probe of the ISM geometry, we find that in general 
the ISM geometries for the majority of the sample appear to 
be somewhere between homogeneous and purely clumpy, im- 
plying that a uniform dust screen geometry does not apply to 
these objects. 

Lastly, we use three different metallicity indicators to mea- 
sure the gas-phase metallicities of the LAEs in our sample. 
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We find that all but one of these galaxies are have lower metal- 
licities than field galaxies at the same stellar mass, with typi- 
cal gas-phase metallicities of ^ 0.2 - 0.4 Zq. Comparing to 
recent results at high-redshift, it appears that galaxies exhibit- 
ing Lya in emission may be systematically more metal poor 
at all redshifts, though more work is needed. 

Studying LAEs spectroscopically allows us to learn about 
their physical properties in more detail and with greater confi- 
dence than with stellar population model fitting. Ideally, these 
low-redshift LAEs would be used as analogs to allow us to 
decipher the physical characteristics of high-redshift LAEs. 
As these objects studied here are more evolved that typical 
narrowband-selected high-redshift LAEs, we can only put 
limits on what we expect to find at high-redshift. However, 
with the numerous multi-object near-infrared spectrographs 
coming online soon, we can expect to make great progress 
in understanding LAEs by directly measuring their physical 
properties at z ~ 2 - 3. 
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